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1 Introduction 



The answer to whether there is any new physics beyond the Standard Model (SM) is probably 
in the affirmative. One of the most promising candidates for new physics (NP), and also most 
studied, is Supersymmetry (SUSY), in both its R-parity conserving (RPC) and R-parity violating 
(RPV) incarnations. Unfortunately for the believers of a Theory of Everything, SUSY introduces 
a plethora of new particles, and even in its most constrained version, a few more arbitrary input 
parameters over and above to that of the SM. Thus, it has become imperative to constrain the 
SUSY parameter space as far as possible from existing data. 

There are two aspects to this practice. First, take the experimental data and find how much 
space we can allow for the SUSY parameters. Second, take the bounds obtained by method 1 and 
see what signals one should observe in present or future experiments (and try to explain if there 
are any apparent anomalies in the present data). Taken together, these two methods form a strong 
tool to observe indirect SUSY (or for that matter, any NP) signals, which is complementary to 
the direct observation of the new particles in high-energy machines like the Large Hadron Collider 
(LHC). 

Since SUSY breaks at a high scale, there is no compelling mechanism to suppress the flavor- 
changing neutral current (FCNC) effects once we take the renormalization group (RG) evolutions 
into account. The problem is most severe for the gravity-mediated SUSY breaking (SUGRA) type 
models. A huge number of models have been proposed to solve this problem; we will not go into 
them here. Rather, we will focus on an equally important area of study, viz., constraining the 
FCNC parameters from experimental data. This shows which models are to stay and which are 
not. 

The problem is worse in the RPV version of SUSY, since there exists a large number of FCNC 
couplings from the very beginning. There are no theoretical limits on these couplings except that 
they better not be nonperturbative even at the scale of the Grand Unified Theories (GUT). The 
only way to constrain the individual couplings and the products of two (or more) of them is from 
experimental data. 

We will just discuss, as a sample process, the mixing of neutral B mesons, and the effects 
of RPC and RPV SUSY on it (more processes will be discussed in a subsequent publication). 
Why this process? The reasons are manyfold: (i) The theoretical part of SM and SUSY are 
both well-known (including higher-order QCD corrections), apart from the uncertainties in some 
of the inputs; (ii) The mixing and the CP-asymmetry data have reached sufficient precision, and 
order-of-magnitude improvements are likely to occur in near future; (iii) The SM amplitude is 
one-loop, so that the RPC SUSY amplitude, which must be a one-loop process, has a fair chance 
of competing, even with high sparticle masses; (iv) The RPV SUSY amplitude is also one-loop 
and competes on the same ground (there may be even trcc-lcvcl amplitudes, but the couplings are 
highly suppressed); (v) There are some decay channels, e.g., B (pK [1], B ^ rj'K [2], B ^ tttt 
[3], which indicate that there may be signals of NP hidden in B decays (although a SM explanation 
is never ruled out). 

Effects of RPC SUSY in — B^ mixing have been exhaustively studied in the literature, and 
constraints were put on the FCNC parameters of different SUSY models [4, 5, 6]. Apart from that, 

SUSY effects on various B decay processes (radiative, leptonic, semileptonic and nonleptonic) have 
been thoroughly investigated, but we are not going to discuss this aspect in the present paper [7]. 
A similar exercise has been performed for RPV SUSY models too [8, 9]. 
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In this paper, we will explore the robustness of some of the FCNC parameters in RFC and 
RPV SUSY as quoted in the literature [5, 6, 9]. The importance of this is threefold: first, this 
will serve as an update of the existing results in the light of new data; second, this will show that 
the bounds can get substantially relaxed once wc take all the SM uncertainties into account (a 
systematic study of the SUSY FCNC parameters, taking this point consistently into account, has 
not been performed as far as we know, though Ref. [6] does it in a sufficiently exhaustive way for 
RFC SUSY only); and third, since the constraints obtained in this paper are the most conservative 
ones, it should tell the experimenters what sort of signal are to be expected at the most. We will 
also show how robust these bounds are if there happens to be a natural cancellation between RFC 
and RFV SUSY effects. 

Here one must note the limitations of the B-related experiments in constraining the NF mod- 
els. Apart from the experimental uncertainties (sin(2/3), CKM elements, branching fractions, etc.), 
there are a number of inherent theoretical uncertainties, most of which stem from the nonpertur- 
bative nature of low-energy QCD. To disentangle signatures of NF, one must be fairly lucky to get 
a sizable deviation from the SM prediction. This is precisely the reason why signatures of RFC 
SUSY will be unobservable in those decay modes which have a tree-level amplitude in the SM. 

We will consider only the gluino-mediated box for the RFC SUSY amplitude. This constrains 
the mixing parameters 5f-^ between the first and the third generations for the down-quark sector. 
Though neutralino and charged-Higgs boson mediated boxes are expected to be small compared 
to the gluino box, an almost equally large contribution comes from the chargino diagram [10, 11]. 
Since the latter diagram constrains the mixing in the up-quark sector (the ^"3 parameters), it will 
not be relevant for our future discussion. The relevant details are to be found in Section 2. For 
RFV SUSY we will consider only one product coupling giving rise to a new mixing amplitude to 
be nonzero at a time. 

We will assume, just for simplicity, that NF affects only the mixing amplitude, but there is 
no NF in the subsequent decay processes. On the one hand this means that the CF asymmetry in 
the channel B J/ipKs measures the phase in the box amplitude. This is easy to implement for 
RFV SUSY by choosing appropriate combinations of nonzero flavor-dependent couplings ^. For 
RFC SUSY, this assumption amounts to neglecting loop-induced SUSY contributions to tree-level 
SM ones, which is a safe assumption. 

The QCD corrections to the box amplitude has been implemented upto Next to Leading Order 
(NLO), both for RFC and RFV SUSY. The relevant anomalous dimension matrix may be found 
in [6]. 

The paper is organized as follows: in Section II, we briefly recall the FCNC phenomena in 
RFC and RFV SUSY, and discuss our numerical inputs in Section III. Section IV deals with the 
constraints coming from RFC SUSY, while Section V does the same job for the RFV version. In 
Section VI we summarize and conclude. 

^However, this is true only for channels like B — > J/tpKs- The channel B — » tt+tt", which measures a, has RPV 
contributions both in mixing and in subsequent decay, and they must be treated together [12]. 
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2 - Mixing in SUSY 



2.1 B« - BO in the SM 



The off-diagonal element in the 2x2 effective Hamiltonian for the neutral B system causes the 
mixing between the gauge eig enstates B^{= hd) and B^{= bd) [13]. The mass difference between 
the two mass eigenstates AM^ is given by 



AMd = 2|Mi2| 



(1) 



where 



Mi2 = 



_ {BO\H,ff\B^) 



2mB 



G 



^{VtdVtlfvBmBfBBBmwSoixt), 



with xt = mf/rn^r, and 



Sq{x) = 



Ax 



llx^ + x^ 



3a;^ In x 



(2) 



(3) 



4(1 -x)2 2(1 

We follow the convention of Ref. [13] for normalization of the meson wave functions. The perturba- 
tive QCD corrections are parametrized by r/^, while the nonpcrturbative corrections are dumped 
in Bb- /b is the B meson decay constant. The sublcading boxes with two charm quarks or one 
charm and one top quark are entirely negligible (the same holds for RPC SUSY; however, due to 
the nonuniversal nature of the relevant couplings, this is not true for the RPV version). 

For B decays to a flavor-blind final state / (e.g., J/tl^Ks) where there are no nonzero CKM 
phases in the decay amplitude, the measured CP asymmetry is proportional to the imaginary part 
of the mixing amplitude. For B^ — B^ box, this is sin(2/3) as arg{Vtd) = — /? (we will implicitly 
assume the Wolfenstein parametrization of the CKM matrix, though the physical observables are 
parametrization invariant). For Bs box, the amplitude is real, so there is no CP violation in the 
SM (to the leading order, i.e., neglecting 0{X^) terms in the CKM matrix, where A = Vus ~ 0.22). 

NP adds up one (or more) new term to M12. Even if it is real, the effective phase /3e// should 
change from (3. Thus, AM^ and Acp{B J/ijjKs), taken together, should constrain both the 
real and the imaginary parts of the NP amplitude. 



Suppose there is one NP amplitude with a weak phase of —2^ so that one can write 



M- 



12 



jMfs^j exp(-2i/3) + \M^2^\ exp(-2i0). 



This immediately gives the effective mixing phase /3e// as 

^ " |Mf2^| cos(2/?) + \Ml^2 \ cos(2(^) 

and the mass difference between the B meson mass eigenstates as 

AMd = 2 f|Mf2^|2 + |Mi^^|2 + 2|Mf2^||Mi^^| cos 2(/3 - , 



1/2 



(4) 



(5) 



(6) 



These are going to be our basic formulae. The only remaining task is to find 



NP 
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Figure 1: Gluino-mediated SUSY contributions to — B'^ mixing. One needs to add the crossed 
diagrams too. The crosses on the squark propagators denote the insertion of the relevant S param- 
eters. All chirality combinations are possible for the quarks, which we do not show explicitly. 

Note that even if ^ = 0, /3e// / /?, which means that even for real NP amphtudes, the 
constraints vary with the choice of (3. The CP asymmetry in S — J/i/^Ks measures /9e//. However, 
with NP contributing to AM^, the B° -W mixing input to the standard CKM fit is lost, and /3 
essentially becomes a free parameter. The same happens for Vtd also. We will discuss these issues 
in Section III. 

2.2 R-Parity conserving SUSY 

In RPC SUSY models there can be two more independent phases <^a and <^b apart from the CKM 
phase, but the electric dipole moment of the neutron constrains them to be small (~ O(10~^-10^^) 
unless the squarks are extremely heavy or there is a fine-tuning betwen them [14]. We take both 
of them to be zero, a choice which can be theoretically motivated, since ^A{(f>B) is the relative 
phase between the common trilinear A-tcrm (bilinear B-term) and the common gaugino mass M 
[15]. Even then one can have new contribution to CP violation coming from SUSY FCNC effects 
[4]. The origin of SUSY FCNC can be easily understood: quark and squark mass matrices are not 
simultaneously diagonalizable. At ~ "mw, radiative corrections induced by up-type (s)quark 
loops are important. These corrections arc typically of the order of log{As/mw) (A5 is the SUSY 
breaking scale) and hence can be large for SUGRA type models. This generates FCNC which 
occurs even in the quark-squark-neutral gaugino vertices, but the flavor structure is controlled 
by the CKM matrix. However, this last feature need not be true in any arbitrary SUSY model, 
particularly those with nonuniversal mass terms. 

One generally works in the basis where the quark fields are eigenstates of the Hamiltonian. 
SUSY FCNC can be incorporated in two ways: (i) Vertex mixing, an approach where the squark 
propagators are flavor and 'chirality' ^ conserving, and the vertices violate them; (ii) Propagator 
mixing, where flavor and 'chirality' are conserved in the vertices but changed in propagators. The 
second approach is more preferred for phcnomenological analysis, since the higher order QCD 
corrections are known better. This is also known as the Mass Insertion Approximation (MIA) 

^Squarks cannot have a chirality, but this is just a loose term to denote the partners of respective chiral fermions. 
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[16, 5]. 

At the weak scale one can write the 6x6 squark mass matrix (say the down type) as 



,2 _/'A^l>LLr^ + AL Ai« 



where the A terms incorporate the FCNC effects. Different FCNC effects are parametrized in 
terms of (5^^ = A^^/m^, where m = ^m\m2, the geometric mean of the masses of the two 
participating squarks, i and j are flavor indices, and A and B are chiral indices. The 5s are 
completely calculable in constrained MSSM, but this is not true in general. Theoretically, for the 
success of perturbative analysis, one expects \5\ < 1. The gluino-mediated box diagrams causing 

— mixing constrain only S^^s of different chiralities in the down-quark sector. (Thus, 
in our notation means {5i-^)ab in the more conventional notation; there is, however, no chance 
of confusion since we do not deal with the 5" parameters.) The same thing can be said for the 
up-quark sector by considering the chargino contributions [10]. 

The AB = 2 effective Hamiltonian can be written in a general form as 

5 3 

'^efT^ = E ^'^^ + E + B.C. (8) 

i=l i=l 

where 

01 = ih^PLd)i{h,.PLd)i, 

02 = {bPRd)l(bPRd)l, 

03 = (bPRd)s{bPRd)s, 

04 = ibPLd)i{bPRd)i, 

05 = {bPLd)s(bPRd)s, (9) 

and Pr(l) = (1 + (~)75)/2- The subscripts 1 and 8 indicate whether the currents are in color- 
singlet or in color-octet combination. The OjS are obtained from corresponding OjS by replacing 
R. 

The Wilson coefficients have been computed at the high scale Ms (chosen to be the arithmetic 
mean of the average squark mass and the gluino mass) by evaluating the diagrams in Fig. 1. We 
quote the results [5]: 

Cl = -R(24xMx)+66h{x)){Sllf, 

C2 = -i?(204x/6(x))(4i)2, 

C3 = i?(36x/6(x))(<5^i)2, 

C4 = -i?{(504x/6(x)-72/6(x)),5ii<5^3«-132/6(x),5iy^i}, 

C5 = -i?{(24x/6(x) + 120/6(x)),5ii5)f^-180/6(x),5iy)f^}. (10) 

Here R = a^/216m| and x = m^/m^. The coefficients qs can be obtained from corresponding CjS 
again with L ^ R. The functions /e and /e are given by 

6(1 + 3x) Inx + - 9x^ - 9x + 17 



feix) 



6{x - If 



~ . . 6.r (1 + ,r) In ,r - .r'^ - 9.r- + 9,r + 1 

/6(^) = ^^^31)5 • (11) 
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Next, one should evolve these coefficients down to the low-energy scale, taken, following Ref. [6], 
to be /X = mb = 4.6 GeV, using the NLO-QCD corrections. The low-scale Wilson coefficients are 

^^{^^) = E E {^r' + ^'^^C,(M5) (12) 

r s 

where i] = as{Ms) / Ois{mt). For the numerical values of a, h and c matrices we refer the reader to 
eq. (10) of Ref. [6]. 

The operators Oi are also to be renormalized at the scale /x. The expectation values of these 
operators between and B^ at the scale /i are given by 

mn)) 

(04(/x)) 

where 

S = 

The S-parameters, whose numerical values are given in Section 3, have been taken from [17]. Note 
that the expectation values are scaled by factor of 2mB over those given in some literature due to 
our different normalization of the meson wavcf unctions. It is trivial to check that both conventions 
yield the same values for physical obscrvables. 

We wish to draw the reader's attention to the fact that with changing x, the interference 
pattern between the SM box and the SUSY box changes. For example, if only ((5j^^)^ is nonzero, 
there is a constructive interference with the SM box if x < 2.1, and a destructive interference 
otherwise. This is just because ci changes sign as we go to higher values of x. Thus, if one chooses 
(^ll)^ to be real, /3e// goes down from /? for low x, and goes up for high x. Near the crossover 
region, the SUSY contribution almost vanishes, so one can in principle have large S parameters. 
We do not analyze these regions since they smell of a fine-tuning, but one should keep this point 
in mind. 



"^b/I-Bi(/x), 



= --5m|/|52(M), 

= ^Sm\flB2,{ii), 

= ^5m|/|S4(/x), 

= ^SmlflB^iii), 



ruB 



(13) 
(14) 



2.3 R-Parity violating SUSY 

R-parity is a global quantum number, defined as (—1)35+^+25^ which is +1 for all particles and — 1 
for all supcrparticles. In the minimal version of supcrsymmetry and some of its variants, R-parity is 
assumed to be conserved ad hoc, which prevents single creation or annihilation of superparticlcs. 
However, models with broken R-parity can be constructed naturally, and such models have a 
number of interesting phenomenological consequences [18, 19]. Some of these R-parity violating 
models can be motivated from an underlying GUT framework [20]. 

It is well known that in order to avoid rapid proton decay one cannot have both lepton number 
and baryon number violating RPV model, and we shall work with a lepton number violating one. 
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Figure 2: R-parity violating contributions to — mixing. Figure (a) corresponds to L4, while 

figure (h) to L2 amplitudes (see text for their meanings). One must add the crossed diagrams, as 
well as the diagrams in (h) where the W is replaced by the charged Higgs or the charged Goldstone. 
The internal slepton can be of any generation, and so can be the internal charge +2/3 quarks, 
generically depicted as u. 



This leads to slepton/sneutrino mediated B decays, and new amplitudes for B^ — B^ mixing 
with charged sleptons and up-type quarks (and maybe W, charged Higgs and charged Goldstone 
bosons) flowing inside the loop (see Fig. 2). Since the current lower bound on the slepton mass is 
weaker than that on squark mass, larger effects within the reach of current round of experiments 
are more probable in this scenario. We start with the superpotential 

Wy = KjkLiQjDl (15) 

where i,j,k = 1,2,3 are quark and lepton generation indices; L and Q are the 5C/(2)-doublet 
lepton and quark superfields and D'^ is the SU (2)-singlet down-type quark superfield respectively. 
Written in terms of component fields, this superpotential generates six terms, plus their hermitian 
conjugates, but for our present purpose the only relevant term is 

CiD-y.^j^e^Ldnui + H.c. (16) 

With such a term, one can have two different kind of boxes, shown in Fig. 2, that contribute to 
^0 _ j^o mixing: first, the one where one has two sleptons flowing inside the loop, along with two 
up-type quarks [21], and secondly, the one where one slepton, one W (or charged Higgs boson 

or Goldstone boson) and two up-type quarks complete the loop [9]. It is obvious that the first 
amplitude is proportional to the product of four A' type couplings, and the second to the product 
of two A' type couplings times Gp. We call them L4 and L2 boxes, respectively, for brevity. 

The effective Hamiltonian for the L4 boxes reads (no sum over k) 

- 32;r2mf ^ J ^^^^ 

where rhi is the slepton mass, and Oi has been defined in eq. (9). Actually, there can be differ- 
ent up-type quarks and different generations of sleptons flowing in the box, but that makes the 
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Hamiltonian proportional to the product of four A' couplings, which we avoid for simplicity. The 
function I{x) is given by 

l-a;2 + 2xloga; 
= ^3^)3 • (18) 

For L2 boxes, there are three different types of amplitudes in the Feynman gauge: involving, along 
with a slepton, a W, a charged Higgs boson, or a charged Goldstone boson. The sum is given by 
[9] 

nL2 = f(l + cot^ (3)xlJixk) + lixk)] O4 (19) 



where Xk = ml /m^ , 



4V2n 

Ukl ""rni 



J[x) = -^i--^)+_i-^+^)^o^- (20) 

and cot/3 = va/vu, the ratio of the vacuum expectation values of the two Higgses that give mass 
to the down- and the up-type quarks respectively (not to be confused with the phase of Vtd)- The 
Hamiltonian is slightly modified if one has different up-type quarks and ui in the box: 

Hl2 = ^^^JI^K^ftK^d [(1 + cot2 f3)x,XpKixk) + Lixk)] O4 (21) 

where we have assumed ruk > nip; if not, the arguments of K and L are to be replaced by Xp. 
The functions are given by 

x-1- logx 

^("^ = ix-ir ' 

L{x) = l-xK{x). (22) 

Note that one can have an imaginary part in the amplitude when the internal quarks are both 
light, but we neglect that effect for our present purpose ^. Also note that in the light of LEP data 
which definitely favours tan /? > 2 — 3 [23] , the Goldstone contributions arc dominant over the 
charged Higgs contributions, which are suppressed by cot^ (3. In deriving the above expressions, 
we have assumed all scalars flowing inside the box to have equal mass. 

In general both the A' type couplings can have a phase, but one of them can be absorbed in 
the definition of the slepton propagator, so it is enough to consider the one remaining phase. 

It is easy to see that the relevant equations (4)- (6) need to be slightly modified to include two 
NP amplitudes when same quarks flow in the loop; this being a trivial exercise, we do not show 
the formulae explicitly. For small values of the corresponding A' couplings, Hl2 dominates H^^, 
but the role may get reversed for large values of the product coupling. Thus, one gets two bands 
in the RPV coupling versus AM^ plane, and our bounds correspond to the outer band. There is 
no such complication when only the L2 amplitude is present. 



3 Numerical Inputs 

The important numerical inputs used in our work is shown in Table 1. A few points are to be 
noted. 

^Even from dimensional arguments, the ratio of the imaginary and the real parts of the mixing amplitude should 
at most be of the order of ml /ml [22], and hence the phase introduced by such an imaginary part can be neglected 
in the analysis. 
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Quantity 


Value 


Remarks 


rriB 


5.2794 GeV 


[24] 


mw 


80.423 GeV 


[24] 


AMd 


0.502 ± 0.007 ps-i 


[25] 


sm(2/3e//) 


0.681-0.784 


[26], at la C.L. 




21.4°-26.0° 


Assumed to be < 90° 


7 


44°-72° 


[27], at la C.L. 




166 GeV 


[27] 




4.23 GeV 


[6] 


mb{mb) 


4.6 GeV 






5.4 MeV 






0.55 ± 0.01 


[27] 




230 ± 28 ± 28 MeV 


[27] 


as{mz) 


0.1172 ±0.002 


[24] 


\Vtd\ X 103 


6.3 - 9.6 


[28], at 95% C.L. 


\Vub\ X 103 


2.49 - 4.55 


[28], at 95% C.L. 



Table 1: Input parameters used for the numerical analysis. 



Unless shown in the table, we have not taken the experimental uncertainty of a quantity into 
account. For example, we have used the central values of the CKM elements, except that of Vtd 
and Vub, in our analysis. They are taken from Ref. [28], extracted without the world average of 
sin(2/3). This is justified since now sin(2/3) itself has a NP contribution, and so should not be used 
to extract the values of the CKM parameters. 

Vtd is determined from AM^; thus, the SM fit for Vtd no longer works when there is one or more 

NP amplitudes to the box. Since the CP asymmetries are controlled by phases not all of which 
are in the CKM matrix, the usual argument of the so-called Universal Unitarity Triangle (UUT) 
[29] does not hold. In essence Vtd becomes a free parameter, only controlled by the unitarity of the 
CKM matrix. To take into account this feature, we have taken the 95% confidence limit (CL) for 
Vtd, extracted without the global average of sin(2/3), for our analysis [28]. The same holds for V^b, 
which contains the phase 7. As pointed out by [6], 7 also becomes a free parameter. We address 
this issue by keeping the range of 7 within la CL quoted by [27], while Vub is varied over its 95% 
CL range. This, we have checked, essentially covers the whole region generated by a narrower 
range of Vub and a wider range of 7 covering to 2tt, with the constraint that the three-generation 
CKM matrix is unitary (i.e., the unitarity triangle should close). The bounds are of the same 
order, but slightly more conservative for the former case. 

The imaginary part of 6s crucially depends on the choice of /3. For our analysis, we have 
varied sin(2/3) between and 1, putting a special emphasis on those values which make (3 f» Peff- 

The most important parameter for the RPC SUSY analysis is the average squark mass, which 
we fix at 500 GeV, effectively neglecting the splitting caused by the SU(2) D-term. The gluino 
mass is varied between 0.3 < m|/m| < 4.0. The bounds more or less scale with the squark mass, 
as can be seen from eqs. (8) and (10). This, however, leaves out some models, e.g., those with an 
extremely light gluino. 
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For RPV SUSY analysis, we take all sleptons to be degenerate at 100 GeV. The bounds on 
the product couplings scale as the square of the slepton mass. The charged Higgs is also assumed 
to be at 100 GeV. As is discussed below, the precise value of the Higgs mass is not a crucial input. 

The value of tan/3 (the SUSY parameter) is fixed at 3, compatible with the lower bound of 
the recent LEP analysis [23]. A glance at eqs. (19) and (21) should convince the reader that the 
bounds are not very sensitive to the exact value of tan P, unless it is small, since the Goldstone 
contributions independent of /? control the show. We have explicitly checked the robustness of the 
bounds with the variation of tan /?. 

The 5-parameters have been taken from [17], and at /i = m^, read 

5i = 0.87(4)+^, B2 = 0.82(3)(4), B3 = 1.02(6)(9), = 1.16(3)1^, B5 = 1.91(4)lf . (23) 
The low-energy Wilson coefficients can be found in [6]. 

4 Results for RPC SUSY 



sin 2/3 


X 




\/\&\ 








0.3 


0.046 


0.021 






1.0 


1.0 


0.099 


0.023 








2.0 


0.27 


0.026 








4.0 




0.031 








0.3 


0.017 


0.0075 


0.0030 


0.0040 


0.732 


1.0 


0.036 


0.0080 


0.0033 


0.0068 




2.0 


0.10 


0.0095 


0.0039 


0.0099 




4.0 


0.090 


0.012 


0.0048 


0.015 




0.3 






0.0078 


0.011 


0.5 


1.0 






0.0088 


0.018 




2.0 






0.010 


0.026 




4.0 


0.235 




0.012 


0.039 



Table 2: Bounds on the S^^ parameters when they are all real. 



The real and imaginary parts of the SUSY amplitude being established through the real and 
imaginary parts of the corresponding Ss, it is easy to find the limits on those real and imaginary 
parts. We perform a scan on the complete range of Vtd and /bVBb, as well as on the SUSY 
phase (p, where generically 6^ = \6'^\ exp(— 2^0), over the range to tt. We demand that AM^ and 
sin(2/3e//) should lie between the values specified in Table 1. The results are obtained for various 
values of /3, even for extreme values like sin(2/3) = or 1. 

Tables 2 and 3 summarize our results. In table 2, we show the bounds on various 6 parameters 
(and their combinations) when they are real, assuming only one to be nonzero at a time. Note 

how the interference pattern changes for Sj^j^ with a: ; it is constructive for small x, and destructive 
for large x. Near x = 2, the crossover occurs, so one can have a large value for 5]^^. A look at 
the respective Wilson coefHcients and hence the interference pattern should help the reader to 
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sin 2/3 


X 




^\Im{5f^Y\ 




^\Im{5f^Y\ 




0.3 


0.046 


0.079 


0.021 


0.035 


1.0 


1.0 


0.10 


0.18 


0.022 


0.039 




2.0 


0.27 


0.27 


0.026 


0.046 




4.0 


0.23 


0.30 


0.031 


0.061 




0.3 


0.077 


0.078 


0.028 


0.030 


0.732 


1.0 


0.16 


0.16 


0.035 


0.038 




2.0 


0.27 


0.28 


0.038 


0.040 




4.0 


0.24 


0.25 


0.051 


0.054 




0.3 


0.081 


0.070 


0.031 


0.027 


0.5 


1.0 


0.17 


0.15 


0.038 


0.032 




2.0 


0.31 


0.23 


0.046 


0.038 




4.0 


0.28 


0.19 


0.053 


0.044 




0.3 


0.084 


0.048 


0.034 


0.021 


0.0 


1.0 


0.18 


0.10 


0.040 


0.026 




2.0 


0.27 


0.28 


0.054 


0.035 




4.0 


0.30 


0.24 


0.060 


0.032 



Table 3: Bounds on the real and the imaginary parts of 5^j^ and 5]^^ 



understand the missing entries. Without tlic CP-asymmctry constraint, nontrivial entries would 
occur everywhere [5]. 

Table 3 shows the real and imaginary parts for (5^^ and (5^^^ for several values of sin(2/3). Due 
to the presence of the nonzero SUSY phase, all entries are nonvanishing. Also, the rise of 5^]^ near 
X = 2 is less pronounced. We do not show the entries for negative values of sin(2/3); they are more 
or less equal with their counterparts for positive sin(2/3). There are a few exceptions, which do 
not change the general result. 

We caution the reader that these bounds are the most conservative ones at the particular 
benchmark points that we have chosen, but by no means signify an impossibility of having larger 
(5s at other points in the SUSY parameter space. Particularly, note that these bounds scale with 
the squark mass. 

5 Results for RPV SUSY 

We explicitly assume that the contributions coming from the RPC SUSY sector vanish if there is 
nonzero RPV interactions. This will be justified post hoc when we discuss a sample case where 
both are present, and there is a possibility of cancellation. 

The strategy is the same as that adopted for the previous subsection. The phase of the RPV 
product coupling is varied between and 2tt while the magnitude of the coupling is assumed to be 
only positive. The range of the scan is kept between the direct product limits, i.e., the limit one 
obtains when one multiplies the individual limits for the two A' components. This limit, as can 
be easily checked [30], is most lenient for the third slepton generation. Only two of the relevant 
product couplings have been bounded from other sources: the product yn^^'i^i has a very stringent 
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A' A' combination 


it;e(A'A') 


7m(A'A') 


Direct bound 


(i31)(i33) 


1.8 X 10"^ 


1.7 X 10-3 


0.202 


(i21)(i23) 


1.2 X 10-3 (0.022) 


1.2(1.4) X 10-3 


0.27 


(ill)(il3) 


2.4(2.6) X 10-3 


2.5(2.6) X 10-3 


2.7 X 10-3 (*) 


(ill)(i23) 


0.016 


0.016 


0.057 


(ill)(i33) 


0.026t 


0.026t 


0.026 


(i21)(il3) 


2.5(3.0) X 10-^ 


2.8 X 10-^ 


0.057 


(i21)(i33) 


0.098 


0.1 


0.23 


(131) (123) 


1.4 X 10" ' 


1.35(1.(3) X 10"-^ 


0.26 



Table 4: Upper limits on the real and the imaginary parts of the relevant A'A' couplings for 
sin(2/3) = 0.732. The numbers in parenthesis show the maximum possible value, for some other 
sin(2/3). All the products, except the one marked with a dagger (f), show Improvements over the 
corresponding bounds obtained from direct product of the bounds of the relevant A's (for i = 3), 
shown in the fourth column. All the numbers in the fourth column are from [30], except the one 
marked with asterisk, which is from [12]. The bounds on the imaginary parts are new. 



bound from tree-level — mixing, of the order of 10"*^ (and hence we do not discuss this 
product further), and the product ^iu^'j^i^ lias been constrained from the measured branching 
ratio and CP asymmetries in the B — tt+tt" channel [12] (marked with an asterisk in Table 4). 
There is no bound on the imaginary parts of the couplings. 



0.002 




-0.001 I 1 1 < 1 < 1 1 1 

-0.002 -0.0015 -0.001 -0.0005 _ _ 0.0005 0.001 0.0015 0.002 

Re Vi3iX,i33 

Figure 3: The real and imaginary parts of A^g^A^gg for various values of sin(2/?). 

Table 4 summarizes our results. Note that almost all the products, apart from those two 
discussed above, and the one marked with a dagger, have been improved, some by orders of 
magnitude. Six among these eight entries have been considered by the authors of [9]. One may 
note that they obtained bounds which have the same orders of magnitude to the ones that we get. 
However, there are several ways in which we have improved upon their calculation, apart from 
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taking the updated data as input. These improvements are: (i) imposition of the CP-asymmetry 
constraint, which was not available at their time (this helps us to obtain the bounds on the 
imaginary parts of the couplings); (ii) incorporation of the NLO QCD corrections, which are, 
anyway, expected to be small — at least they should not change the numbers by, say, a factor of 
two; (iii) scan over the full range of the input parameters, as we have already discussed; and (iv) 
consideration of the SM contribution, and the possibilities of interference between SM, L2 and L4 
amplitudes. Also note that they have computed the bounds for tan/3 = 1. 

Our bounds are more or less insensitive to the precise value of sin(2/3) (the UT angle); the 
exceptions arc shown in the table. Another feature is that real and imaginary parts have almost 
the same bounds. Figure 3 highlights this nature in detail for the product coupling X'i^iX'i^^- 

Let us try to understand this figure. There arc three main regions, the first one (like a smeared 
cross) encompassing the origin, the second one, divided into four almost symmetrical fragments, 
around |^?e(A^3iA^33)| = 5 x 10~^, and a third fragmented region about |A^3iA^33| = 15 x 10~^. It 
is the outermost third region which gives the bound. It is easy to explain the origin of these three 
regions, and sheds light on the role of sin(2/3) in determining the parameter space. 

The first one is governed by the SM, and the difference between the experiment and the theory 
is filled up by RPV. Note that this region has points only for sin(2/3) = 0.732; thus, SM is allowed 
without any NP. There is a satellite region, for sin(2/3) = 0.5, where SM is not allowed (from the 
Acp constraint), and RPV fills in to generate the necessary CP asymmetry. The second region 
includes points for extreme values of sin(2/3): or 1. SM is not allowed, and one needs a greater 
role from RPV to obtain the observed CP asymmetry. Though both L2 and L4 boxes are allowed 
for this RPV coupling, the L2 box dominates here. However, they come with opposite sign, so 
there is a region where they interfere destructively (particularly with the increase of the coupling) , 
and RPV contribution may even go to zero, leaving only SM. This generates the third region and 
explains why one has points even for sin(2/3) = 0.732. 

It is a possibility that both RPC and RPV SUSY are present, however pathological that may 
seem. Even in this extreme case the bounds are never changed by orders of magnitude, unless 
there is a very precise fine-tuning. For example, with only and A^2i'^i33 present (this particular 

combination is chosen since both of them have comparable upper bounds), the bounds are relaxed 

to ^\Rei5iir\ = 0.22, ^\Im{5lir\ = 0.26, Re{X',^,\',^,) = 0.12, Im,{X[^,X[^^) = 0.11, for x = 1 
and sin(2/3) = 0.732. Thus we have reasons to be confident about these bounds. This is more so 
if the limits have different orders of magnitude. 



A' A' combination 


Decay channels 


A' A' combination 


Decay channels 


(131) (133) 
(ill)(i23) 
(i21)(il3) 


b diiii, b dVlUi 
b cud, b sdd 
b — > ucd, b — > dsd 


(121) (123) 
(ill)(il3) 


b ccd, b ssd 
b diiii, b diTiUi 

b uud, b ddd 
b — >■ diiii, b dlT^Vi 



Table 5: The possible decay channels of the B meson driven by the different RPV product cou- 
plings. The final state mesons are not shown explicitly. For the semileptonic decays, the outgoing 
leptons must be of the same generation, denoted by i. 

Most of these product couplings contribute to various B-decay channels, both nonleptonic and 
semileptonic. They are enlisted in Table 5. It is easy to check which mesons and leptons come out 
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in the final stage. The present e~^e B factories, as well as the future hadronic machines, should 
put tighter constraints on these RPV product couplings. 

6 Summary and Conclusions 

In this paper we have enlisted the constraints on the real and the imaginary parts of the FCNC 
parameters of both R-parity conserving and R-parity violating SUSY, coming from — 
mixing. For RFC SUSY, these are the conventional parameters of different chiralities. The 
same analysis was performed by [6]; our results differ slightly from theirs due to two reasons. We 
have performed a scan over all SM quantities, including Vta and /^Bb, and our range of scan for 

7 is different from theirs. Weaker constraints on these S parameters can also be derived from the 
radiative decay b d^f. 

For the RPV SUSY scenario, the FCNC parameters are the A' type lepton-number violating 
couplings. One needs a product of two such couplings to generate B^ — B^ oscillation. There can 
be several such products, depending on the choice of the quarks and sleptons flowing inside the 
box, whose real and imaginary parts have been bounded from the experimental data. The bounds 
on the real parts update the work of [9] while the bounds on the imaginary parts are derived for 
the first time in this paper. 

If we ever find a signal for NP in B factories, how can we be certain that it indeed comes 
from SUSY? There are three steps to ascertain that. First, sort out those channels which show 
an abnormality. Second, Try to find the model which can explain these anomalies. And third, 
check whether there arc other channels where one may expect to see an anomaly, and whether 
the anomaly may be present in the data. If there is a prospective channel, one should look for 
it. Confirmation of the nature of the NP is never possible without the study of such correlated 
signals. Such correlated signals may even be the direct production of new particles, e.g., in the 
Large Hadron Collider. A possible discriminating signal between RFC and RPV SUSY is the 
fact that the RPV version may be highly fiavor-specific, and so one would expect the absence of 
anomaly in such channels which may be affected in a more flavor-blind model such as the RFC 
SUSY. 

For the help of our experimentalist colleagues, we enlist the possible decay modes of the B 
meson which are driven by the RPV product couplings discussed here. A careful study of the 
possible channels in present and future B factories should be able to put tighter constraints on 
the parameter space. Particularly, the proposed higher luminosity e^e^ B factories should make 
both sin(2/3) and AM,^ precision observables, and the bounds are expected to be improved by at 
least one order of magnitude, if we do not discover SUSY by that time. 
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